Role of disorder in the size-scaling of material strength 
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We study the sample size dependence of the strength of disordered materials with a flaw, by 
numerical simulations of lattice models for fracture. We find a crossover between a regime con- 
trolled by the fluctuations due to disorder and another controlled by stress-concentrations, ruled 
by continuum fracture mechanics. The results are formulated in terms of a scaling law involving a 
statistical fracture process zone. Its existence and scaling properties are only revealed by sampling 
over many configurations of the disorder. The scaling law is in good agreement with experimental 
results obtained from notched paper samples. 
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The fracture strength of materials depends on various 
characteristic length-scales of the specimen, and repre- 
sents a fundamental open problem of science and engi- 
neering. Probably the oldest scientific study of this issue 
was performed by Leonardo da Vinci, who measured the 
carrying-capacity of metal wires of varying length 
The simple observation was that the longer the wire, the 
less weight it could sustain. The reason for this behavior 
is rooted in the structural disorder present in the mate- 
rial: the strength is dominated by the weakest part of the 
sample and its distribution is related to extreme value 
statistics Longer wires are more likely to contain 
weak parts and are thus bound to fail earlier on average. 
While the physical mechanism behind this extreme- value 
based statistical size effect is clear, obtaining mathemat- 
ical laws for it is still a formidable task. In quasi-brittle 
materials, such as concrete and many other composites, 
this issue is particularly important and it is complicated 
by the significant damage accumulation preceding sample 
failure. 

The most important setting to study size effects in- 
volves a specimen containing a pre-existing flaw, a notch. 
Failure in this case is determined by the competition be- 
tween deterministic effects, due to the stress enhance- 
ment created by notch, and the stochasticity coming from 
the disorder [3[. For a sufficiently large notch, stress en- 
hancement around the crack tip dominates the process. 
It is then customary to consider disorder as a small per- 
turbation, by defining a Fracture Process Zone (FPZ) 
around the crack tip, where all the damage accumulation 
is confined. For quasi-brittle materials, however, the size 
of the FPZ may not be negligible when compared to the 
system size. Conversely, for small notches failure is influ- 
enced by statistical effects, and may for instance initiate 
far from the pre-existing notch due to nucleated microc- 



racks. Several formulations to account for the size effects 
have been proposed in the literature 0, @, H, 0, H, H|. 
These approaches are mainly based on the ad hoc exten- 
sions of linear elastic fracture mechanics (LEFM) . This is 
a well established framework to understand cracks in ho- 
mogeneous media but encounters fundamental problems 
when disorder is strong and homogenization methods are 
not applicable. 

In LEFM the stability of a flaw against failure is given 
by the Griffith's energy criterion for the critical stress, 
a c ~ K c / y/ao, where clq is the linear size of the crack 
and the critical stress intensity factor K c ~ y/EG c is 
a function of the fracture toughness G c and the elastic 
modulus E A scaling law for the size-effect has 

been proposed by Bazant for quasi-brittle materials jij. 
One generalizes the Griffith expression by postulating an 
additional length-scale £ due to the presence of a FPZ 
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Equation JT]) tries to incorporate two natural, impor- 
tant effects: First, in the large notch limit £/ao <C 
one should recover an expression that follows the LEFM 
scaling, in which the strength is inversely proportional to 
1/y/aQ. Second, for a vanishing external flaw size ao — > 0, 
the average strength should still remain finite. 

Here we investigate the role of the disorder in the fail- 
ure of notched quasi-brittle specimens, providing a micro- 
scopic justification and establishing the limits of validity 
of Eq. (JTJ. We study the size scaling of strength by 
the extensive numerical simulations, a difficult task due 
to the different length scales involved and to the need 
of significant statistical averaging. We vary the disor- 
der, which we model as a locally varying random failure 
threshold, and show that it plays a crucial role in deter- 
mining the size effect, influencing the fracture toughness 



K c . Furthermore, the lengthscale £ naturally emerges 
from the simulations and can be shown to be directly 
related to the FPZ size. Finally, for notch sizes smaller 
than a critical length a c , we observe a cross-over to the in- 
herent, sample-size dependent strength of the unnotched 
sample. We present a scaling formula that incorporates 
all these effects and confirm its validity by comparing the 
simulations to experiments on notched paper samples. 

To simulate a disordered elastic solid, we consider the 
simplest case where disorder and LEFM-like stress en- 
hancements can be incorporated. We perform extensive 
simulations of the random fuse model (RFM) [ll|, 
The RFM represents a quasi-brittle failure process by 
an electrical analogue composed of a network of fuses. 
We consider a triangular lattice of linear size L with a 
central notch of length ao. The fuses have unit conduc- 
tance (which would correspond to E = 1 in the elastic 
system) and random breaking thresholds i c . These repre- 
sent a locally varying fracture toughness/strength. The 
i c lie between and 1, with a cumulative distribution 

P{i c ) — ic , where D represents a quantitative measure 
of disorder. The larger D is, the stronger the disorder. In 
the simulation, the burning of a fuse occurs irreversibly, 
whenever the electrical current in the fuse exceeds its 
threshold i c . An external current is increased applying 
a voltage difference between the top and the bottom lat- 
tice bus bars and applying periodic boundary conditions 
along the other direction. Damage accumulates until a 
connected fracture path disconnects the network and one 
can define the strength a c as the total peak current di- 
vided by the length of the bus bar. We also perform nu- 
merical simulations for the random spring model (RSM) 
O , similar to the mesoscale models used routinely for 
The RSM is similar to the RFM but the 
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concrete 

fuses are replaced by elastic springs that break when their 
elongation reaches a random threshold. While the RSM 
represents more faithfully the elastic continuum, the sta- 
tistical properties of the fracture process are analogous 
to those of the RFM 

The inset of Fig. Q] reports the strength, averaged over 
different configurations, with varying ao, D, and L. The 
most instructive way of plotting is to consider the in- 
verted square strength, 1/cr^. Assuming Eq. (p}, it is 
clear that 1/ct^ should become a linear function of ao for 
large enough notches. Plotting the data in such a manner 
in Fig. [1] reveals four interesting features: (i) for ao ^> 1, 
the scaling of Eq. {U is recovered asymptotically; (ii) the 
linear part of the data when extrapolated towards ao = 
reveals a disorder-dependent intercept £(-D), that should 
be related to the size of the FPZ; (Hi) the slope of the 
linear part of the data (1/K^(D)) is disorder-dependent, 
which implies a disorder-dependent fracture toughness 
G C (D); and finally, (iv) a careful observation (see Fig. [2|) 
reveals that for small ao less than a critical crack size a c , 
the strength scaling crosses-over from a stress concen- 
tration dominated LEFM scaling (Eq. [TJ to a disorder 
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FIG. 1: Inset: The strength <r c as a function of the notch 
size oo for several disorders D and system sizes L. The main 
figure shows the scaling plot of the data according to Eq. |T} . 



dominated scaling. In particular, the data presented in 
the inset of Fig. [5] indicates that the strength of the un- 
notched system (for ao = 0) is smaller than the LEFM 
limit K c j\f^ (where is the slope and £ is the in- 

tercept of the lines in Fig. [1] based on Eq. ((T|)). Hence, 
in order to extrapolate the typical sample-size depen- 
dent strength from notched experiments it is necessary 
to know a c . As can be seen, these features are also ex- 
hibited by RSM results. 

The cross-over at a critical crack size a c marks the im- 
portant role of structural disorder or internal damage on 
the size-effect, and relates to the size-effect in unnotched 
samples. For the RFM, the size-effect without a notch 
has been shown by simulations and theoretical arguments 
to have a logarithmic dependence on the linear system 
size, L T3, if 



For many engineering materials, one re- 
sorts to the Weibull theory [l(| as an empirical starting 
point. 

This cross-over from LEFM dominated strength scal- 
ing to disorder dominated strength scaling is illustrated 
in Fig. [3j where we compare the simulation data to ex- 
perimental results on paper samples with varying cen- 
ter notch sizes. The paper data is from Ref. 17]. The 
two data sets presented are for strips of L = 15cm 
cut from laboratory-made handsheets, of fine paper-type, 
with center defects of nominal sizes of ao =0.5, 1, 1.5, 2, 
and 2.5 mm. Tensile tests were performed on 100 sam- 
ples, in order to average the results. In Fig. [3J we observe 
clearly the presence of a c , below which the notches have a 
negligible role on a c , and for crack sizes ao larger than a c 
fracture is ruled by LEFM. We have also fitted the sim- 
ulation data to these experiments, by considering simi- 
lar ao/L ratios and varying the D to obtain reasonable 
agreement with the experimental data. 

The observations from Figs. [1] and can be summa- 
rized into a single scaling theory by noticing that there 
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FIG. 2: A close-up of the strength data for D — 0.6 and 
various ao and L. The inset compares fracture strength in 
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FIG. 3: Comparison of numerical results with the experimen- 
tal data for two kinds of paper. Figure presents size effect for 
small notch sizes using a RFM simulation of a system of size 
L = 128 and a disorder of D = 0.6. 



must be a scale a c above which the LEFM holds and a c 
follows Eq. (P) . For ao <C a c the strength scaling deviates 
significantly from Eq. ([1]) and saturates to a value that 
depends only on disorder and the sample size, a(L,D). 
This is the strength of the unnotched system of size L and 
disorder D. The cross-over from LEFM scaling (Eq. ([1])) 
occurs at a notch size a c which can be obtained from 
Eq. (HJ as l/a(L,D) 2 ~ (a c + £)/K 2 . lt is possible to de- 
scribe this crossover by the following scaling form, valid 
for all ao such that 

% =Z + a f(a c /a ) (2) 



where the scaling function f(y) fulfills the limits 

The length-scale a c ~ (K c (D)/a(L, D)) 2 - £(£>) cor- 
responds to a cross-over scale below which material 
strength is governed by disorder strength and system size. 
When the notch size exceeds this cross-over length scale 
(ao > a c ), fracture is governed by LEFM. At fixed L, 
stronger disorder will increase a c since the strength of 
unnotched samples decays faster than the decrease in K c . 
This is understandable since the disorder masks more ef- 
ficiently the stress concentration due to the notch. Like- 
wise, at fixed D, a c increases with increasing L. 




FIG. 4: The damage (fraction of failed elements or fuses) at 
maximum stress a c for L = 128 and D — 0.6. The black mark- 
ers illustrate the broken fuses in one single, randomly chosen 
sample whereas the color background represents the damage 
averaged over N = 2000 number of samples. The colour code 
indicates the intensity of damage. A damage cloud that in- 
creases in density is clearly visible close to the initial notch. 
The figure demonstrates the screening of damage due to free 
crack faces. It is clear that in the disordered samples, FPZ 
is a statistical zone, visible only when damage profiles are 
averaged over many samples. 

Simulations of the RFM allow us to access the damage 
evolution prior to failure and can thus be used to visual- 
ize the development of the FPZ. As shown in Fig. [4] for a 
single realization of the disorder, at maximum stress we 
only see diffuse damage, without apparent localization 
so that the FPZ can not be observed. When we average 
the damage over different configurations, however, a clear 
FPZ emerges in front of the crack tip (see Fig. 3}. Hence 
the FPZ should be considered a statistical concept, vis- 
ible only when averaging over disorder, while the effect 
of the FPZ is nevertheless seen in the size-effect (£). To 
measure the FPZ size, we consider a projection of the 
average damage along the crack direction and obtain a 
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FIG. 5: a): Damage profiles along the crack axis for var- 
ious disorders D. Damage profiles follow an exponential 
decay on a uniform damage background, i.e., d(x) — A + 
B exp (— x/^fpz), where A and B are constants and x is the 
distance from the crack tip along the crack axis, b): £ vs. 
£,fpz for various disorders D. 



by LEFM and is furthermore seen in experiments. All 
the regimes are summarized in a generalized scaling ex- 
pression (Eq. @) for the strength of disordered media. 
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profile that is decaying exponentially towards a homoge- 
neous background value: d(x) — A + Bexp (— x/^fpz) 
(see Fig. [5^,). We have analyzed the data for different 
values of D and checked that the profiles do not depend 
on ao as long as this is not too far from a c . The LEFM 
stress intensity factor would indicate a l/y/r -like diver- 
gence of the stress at the crack tip. It is evident that 
the observed exponential shape of the damage profile d 
is in contrast to a 1/y/r -like decay and should be nat- 
urally interpreted as a screening of the crack tip caused 
by the disorder. In Figure 03, we plot the deduced frac- 
ture process size £,fpz{D) against the intrinsic scale £ 
that one obtains from the fits of the strength data to 
the Eq. (JT|) and which also is an important part of the 
scaling theory presented in Eq. ©. It can be seen that 
these are linearly proportional indicating that £ is indeed 
a direct measure of the FPZ size. Notice that an expo- 
nential damage zone has been indeed measured in paper 
samples and the corresponding lengthscale was compared 
with the one obtained from Eq. |T]) |l8j |. 

In conclusion, we have resorted to simulations of sta- 
tistical fracture models to analyze the problem of the 
size-effect in the failure of materials. For large notches, 
the simulations recover the expected scaling of LEFM Q 
and allow to relate the effective FPZ size £ to the actual 
average damage profiles. As the notch size is decreased 
we observe a crossover at a novel scale a c to a disorder- 
dominated size-dependent regime that is not described 
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